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EXECUTIVE  SUMMARY 

One  of  the  architectures  being  Investigated  for  the  future  Defense 
Communications  System  an  Integrated  (voice/data),  common  user, 
communication  system.  There  are  several  technical  approaches  to  im- 
plementing such  a system.  These  approaches  involve  different  mixes  on 
such  network  level  issues  as  : analog  or  digital  transmission;  and  circuit, 
packet,  or  hybrid  switching  technology.  The  network  level  issues  cannot 
be  satisfactorily  addressed  until  accurate  performance  models  are  avail- 
able for  the  individual  switching  subsystems  of  the  network. 

Accurate  performance  models  for  hybrid  switching  technology  subsystems 
turn  out  to  be  exceptionally  difficult  (perhaps  impossible)  to  construct 
on  theoretical  grounds  alone.  Therefore,  simulation  models  are  required. 
This  report  addresses  the  construction  and  use  of  a simulation  model  for 
a particular  type  of  hybrid  switch  known  as  the  SENET  switch. 
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I.  INTRODUCTION 


This  note  discusses  possible  types  of  the  Slotted  Envelope  NETwork 
(SENET)  transmission  concept  and  describes  a simulator  for  a subclass 
of  these  types.  By  using  this  simulator  the  performance  of  SENET 
systems  may  be  determined  for  various  traffic  loadings. 

The  Defense  Communications  Engineering  Center  (DCEC)  is  exploring 
various  architectures  for  the  future  Defense  Communications  System 
(DCS).  Among  the  architectures  being  investigated  is  the  Slotted 
Envelope  NETwork  (SENET)  proposed  by  Drs.  Coviello  and  Vena  [1].  The 
SENET  transmission  concept  integrates  both  circuit-switched  and 
packet-switched  traffic  into  a single  integrated  network.  In  the 
SENET  concept  transmission  time  is  first  broken  into  equal  time 
Intervals  called  frames.  Within  these  frames  circuit-switched  traffic 
is  assigned  to  the  first  region  of  the  frame  and  packet-switched 
traffic  to  the  second.  The  next  section  of  this  note  discusses  the 
various  ways  traffic  may  be  allocated  to  these  regions;  each  of  the 
allocation  techniques  determines  a SENET  subtype. 

The  note  then  describes  a simulator  which  simulates  a class  of 
these  subtypes.  The  class  of  SENET  subtypes  includes  those  which  have 
either  a free  or  upper  constrained  region  for  the  circuit-switched 
traffic;  which  allocate  circuit-switched  traffic  to  a fixed  slot  size; 
and  which,  when  inactivity  Is  detected  within  a circuit-switched 
channel,  compress  the  region.  The  circuit-switched  traffic  is 
assumed  to  be  all  voice  conversations.  Silences  in  conversational 


speech  are  detected  as  inactivity.  A three  state  model  simulates 
conversational  speech  consisting  of  periods  of  talking  followed  by 
either  short  or  long  periods  of  silence. 

The  report  is  divided  into  four  main  sections  and  an  appendix. 
Section  II  discusses  the  subtypes  of  SENET.  The  actual  simulator 
code  is  described  in  the  next  section,  and  the  final  section  describes 
the  Job  Control  Language  (JCL)  cards  required  to  run  the  simulator, 
statistics,  and  plotting  programs.  An  example  is  also  given  of  the 
JCL  needed  to  translate,  compile,  link,  and  load  another  model. 
Finally,  the  appendix  describes  the  development  of  the  speech  activity 
model  used  in  the  simulator. 

The  work  reported  in  this  Technical  Note  is  part  of  the  author's 
dissertation  for  the  degree  of  PhD  from  the  Air  Force  Institute  of 
Technology. 


2 


II.  SENET  SUBTYPES 


The  Slotted  Envelope  Network  (SENET)  [1]  is  a technique  which 
integrates  real-time  data  (digitized  voice,  facsimile,  etc)  and 
packet-switched  data  into  a single  digital  communications  network. 
Basically,  the  system  is  a dynamically  allocated,  synchronous,  time 
division  multiplex.  Transmission  time  on  links  is  first  partitioned 
into  equal  size  periods  called  frames. 

Each  frame  is  further  divided  into  a start  of  frame  and  two 
regions.  The  start  of  frame  is  a series  of  bits  used  to  maintain  frame 
integrity.  The  remaining  two  regions  contain  the  network  traffic.  The 
first  region  contains  the  real-time  traffic  and  the  second  region 
contains  packet-switched  traffic.  The  method  of  allocating  data 
within  each  region  varies. 

Region  I is  subdivided  into  "slots"  and  each  real-time  channel 
is  assigned  a slot.  The  SENET  concept  does  not  specify  that  these 
slots  be  a fixed  or  a variable  size.  Nor  does  the  SENET  concept 
specify  that  a real-time  channel  will  occupy  a single  slot  or  many. 

The  assignment  within  region  II  also  allows  some  room  for  interpreta- 
tion. 

The  second  region  is  used  for  packet-switched  data.  If  there  is 
no  data  in  this  category,  "idle"  characters  will  be  sent  on  the 
transmission  link  until  the  next  start  of  frame.  A data  packet 
usually  contains  a length  count  as  part  of  its  header;  therefore  slots 
are  not  needed  in  this  region.  A problem  arises, however,  when  the 
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TABLE  I.  SUBTYPES  OF  SENET 


I.  BOUNDARY  MOVEMENT 

A.  Fixed 

B.  Movement 

1.  Free 

2.  Constrained 

a.  Time 

(1)  Fixed 

(2)  Varying 

b.  Limits 

(1)  Upper 

(2)  Lower 

(3)  Both 

II.  REGION  I ALLOCATION 

A.  Slot  Size 

1.  Fixed 

2.  Variable 

3.  Incremental 

B.  User  Activity  Detection 

1.  None 

2.  Fill  slots  with  region  II  data 

3.  Compress  region  I 

4.  Use  for  other  region  I data 
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space  remaining  in  the  region  is  too  small  to  contain  a full  packet. 
Two  choices  are  available: 

• Hold  packet  until  start  of  region  II  in  next  frame. 

• Transmit  part  of  the  packet  this  frame  and  the  remainder  at  the 
start  of  the  region  II  in  the  next  frame. 

The  second  method  is  preferred. 

In  addition  to  the  allocation  strategies  within  the  regions  the 
allocation  of  the  two  regions  within  the  frame  may  vary.  If  we  assume 
a "boundary"  between  the  two  regions  then  the  method  of  partitioning 
the  frame  into  the  two  regions  dictates  how  this  boundary  moves.  In 
general  the  boundary  may  be  fixed,  freely  movable,  or  constrained. 

In  the  fixed  boundary  case,  a fixed  portion  of  the  frame  is 
allocated  to  region  I and  a fixed  portion  to  region  II.  In  this  case 
unused  bits  within  the  region  I must  be  flagged  as  idle.  The  portion 
of  the  frame  containing  these  idle  bits  is  effectively  lost. 

In  the  free  boundary  case,  channels  of  the  region  I type  are 
allocated  slots,  as  needed,  up  to  the  entire  frame.  When  no  real-time 
channels  are  present,  there  is  no  region  I.  Therefore,  an  unused 
portion  of  region  I is  not  lost  as  in  the  fixed  boundary  case. 

The  constrained  boundary  falls  between  the  fixed  and  free  cases. 
There  are  various  types  of  constrained  boundaries: 

• Upper  limit, 

• Lower  limit, 

• Range  limits. 


J 
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As  their  names  imply,  the  constraints  may  be  an  upper  and/or  lower 
bound  on  the  portion  of  the  frame  allocated  to  one  of  the  regions.  In 
addition,  these  constraints  may  be  fixed  in  value  or  allowed  to  vary 
with  time  (as,  for  example,  dependent  on  overall  network  loading). 

To  summarize,  SENET  is  a time  division  multiplex  for  real  time  and 
packet-switched  data.  It  handles  each  type  of  traffic  in  one  of  two 
regions  within  its  frame.  Traffic  within  region  I is  assigned  space 
on  a slot  basis  and  within  region  II  on  a serial  basis.  The  boundary 
between  the  regions  may  be  fixed  or  movable.  If  it  is  movable  then  it 
may  be  either  free  or  constrained.  For  the  constrained  boundary  the 
limits  may  be  fixed  or  time  varying. 

Another  dimension  of  classification  is  the  allocation  within 
region  I.  Users  are  allocated  channels  within  region  I for  as  long  as 
the  channel  is  needed.  When  the  user  leaves  the  system,  its  alloca- 
tion is  placed  in  a free  status.  Channels  are  allocated  physical 
space  in  region  I on  a fixed,  variable,  or  incremental  basis.  In  fixed 
space  allocation,  slots  within  the  region  I are  of  a predetermined 
size.  Not  all  slots  need  be  exactly  the  same  size.  When  a user 
requests  a channel,  the  region  I free  slots  are  scanned  until  a slot 
of  the  required  size  is  found.  The  variable  slot  size  is  the  opposite 
extreme.  Slot  sizes  are  not  predetermined;  rather,  they  are  dynamic- 
ally determined  based  on  the  user  requirements.  This  method  becomes 
unwieldy  in  practice,  especially  for  users  requesting  channels  with 


rates  having  fractional  bits  per  frame.  The  Incremental  slot  size  Is 
a compromise  between  the  extremes.  Slots  may  be  of  any  multiple  of  a 
given  minimum  size. 

Another  factor  in  the  allocation  techniques  for  region  I is 
activity  detection.  Some  users  need  a guaranteed  time  slot,  but  may 
not  need  it  all  the  time.  This  is  the  case  for  voice  circuits  with 
silence  detection  and  for  certain  command  and  control  circuits  which 
transmit  a short  "Are  you  alive?"  query  once  every  few  seconds. 

Assuming  that  lack  of  activity  is  detected,  the  slot  may  be  assigned 
for  other  purposes  during  periods  of  "silence".  The  first  alternative 
is  to  "stuff"  the  slot  with  bits  which  normally  would  be  transmitted 
during  the  following  region  II.  (This  has  been  loosely  called  Time 
Assignment  Digital  Interpolation,  TADI,  when  referring  to  digital 
voice).  Another  alternative  is  to  compress  the  region  I portion  of  the 
frame  when  there  is  silence,  and  correspondingly  expand  region  II.  A 
third  alternative  is  to  overbook  the  region  I,  as  Time  Assignment 
Speech  Interpolation  (TASI)  is  used  on  submarine  cables. 

The  allocation  methods  for  region  I traffic  also  determine  sub- 
classes of  the  SENET  concept.  Therefore,  the  method  of  assigning 
slots  to  channels  in  the  region  I portion  of  the  SENET  frame  may  be 
based  on  a fixed,  varying,  or  incremental  slot  size,  and  the  activity 
on  the  channels  may  or  may  not  be  detected.  If  activity  detection  is 
used,  the  detected  unused  capacity  may  be  used  for  data,  compressed 
from  the  region,  or  used  by  other  region  I channels. 
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III.  SIMULATION  MODEL 


The  SENET  simulation  program  Is  detailed  in  this  section.  In 
order  to  help  the  reader  understand  the  model,  a brief  functional 
description  of  the  simulator  is  given.  The  Simulation  Oriented 
Language  ( SOL-370)  used  to  code  the  model  is  then  described.  The  next  • 
part  of  the  discussion  describes  the  values  input  to  the  simulator  and 
the  files  used.  Following  this,  the  actual  simulation  code  is  described. 
1.  SENET  SIMULATION 

The  SENET  concept  integrates  circuit-switched  and  packet-switched 
traffic  onto  a single  transmission  link.  The  simulator  models  a 
single  link.  The  traffic  which  loads  the  link  has  the  following 
characteristics: 

• All  circuit-switched  traffic  is  digitized  voice. 

• Voice  activity  may  be  detected. 

• Voice  calls  have  a Poisson  arrival  with  an  interarrival 
mean  time  of  VIN  ms  (milliseconds). 

• Voice  calls  have  an  exponential  hold  time  of  VHOLD  seconds. 

• Packet-switched  messages  have  a Poisson  arrival  with  an 
interarrival  mean  time  of  DIN  ms. 

• Packet-switched  messages  have  an  exponential  length  with 
mean  DLEN  bits. 

The  simulator  attempts  to  be  a good  representation  of  the  SENET 
link  allocation,  but  many  overhead  control  functions  that  would  exist 
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in  an  actual  switch  are  ignored.  No  overhead  bits  for  the  start 
of  frame,  call  slot  reservation,  or  packet  headers  are  included  in 
the  simulated  SENET  frame.  Calls  entering  the  system  have  no  setup  or 
teardown  times.  The  number  of  bits  in  each  message  is  added  to  the 
queue,  but  no  attempt  is  made  to  identify  individual  messages  or 
packets.  All  bits  transmitted  in  the  packet-switched  data  region  of 
the  frame  are  removed  from  the  data  queue  during  the  next  frame. 

The  simulator  does  try  to  emulate  realistically  data  arrival 
and  frame  allocation.  Voice  calls  and  messages  arrive  randomly.  If 
there  is  room  for  its  allocation,  the  voice  call  is  immediately  added 
to  the  number  of  calls  in  progress,  NVOICE;  its  ending  time,  TEND,  is 
determined;  and  its  voice  activity  is  sensed.  As  soon  as  a message 
arrives  its  length  in  bits  is  determined  and  added  to  the  data  queue, 

QUEUE. 

The  simulator  assumes  the  frame  has  a time  period  of  TFRAME  ms 
and  a size  of  MSFRAME  bits.  The  circuit-switched  region  has  an  upper 
constraint  of  MVFRAME  bits  which  must  be  large  enough  to  allocate  up 
to  MNVOICE  calls  with  a fixed  slot  size  of  VSLOT  bits  each.  If  a call 
arrives  that  would  cause  the  number  of  "i n-progress"  calls,  NOVICE,  to 
oe  more  than  MNVOICE  calls,  the  arriving  call  is  lost  and  the  lost 

calls  counter,  LVOICE,  is  incremented. 

The  frame  is  allocated  by  the  simulator  in  the  same  time  sequence 
as  in  the  SENET  concept.  The  time  of  the  next  start  of  frame, 

NEXT  SOF,  starts  the  frame  allocation.  The  speech  activity  during  the 
previous  frame  of  each  in-progress  call  is  determined.  The  allocation 
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for  the  call  in  the  circuit-switched  region,  VFRAME,  of  the  SENET 
frame  is  changed  if  needed  Any  calls  that  ended  during  the  previous 
frame  are  removed  from  the  number  of  calls  in  progress,  NVOICE,  and 
any  present  allocation  in  the  circuit-switched  region,  VFRAME,  is 
removed. 

Data  bits  transmitted  during  the  previous  frame,  NDATA,  are 
removed  from  the  data  queue,  QUEUE.  Then  the  number  of  bits  which  may 
be  transmitted  this  frame  is  calculated  and  added  to  the  region  I 
allocation,  VFRAME,  to  form  the  total  SENET  frame  allocation,  SFRAME. 

The  above  process  repeats  until  the  simulation  time,  TIME, 
exceeds  the  simulation  time,  SIMT,  specified  by  the  user. 

2.  THE  SOL-370  LANGUAGE 

The  SENET  simulation  model  is  written  in  SOL-370  (release  1-79). 
This  simulation  language  was  chosen  because  of  its  simplicity  and 
and  availability.  This  language  is  an  extension  of  the  language 
described  by  Guffee  and  Ulfers  [2]. 

SOL  was  originally  developed  by  D.  E.  Knuth  and  J.  L.McNeley  [3] 
and  initially  implemented  on  a Burroughs  B-5500.  Because  the  original 
implementation  of  SOL  was  written  in  ALGOL,  which  is  not  supported  on 
the  IBM  370,  DCA  implemented  a dialect  written  in  PL/I.  To  provide 
the  additional  capabilities  needed  for  this  model,  I expanded  the 
language. 

Within  SOL  the  basic  entity  which  "flows"  in  the  simulation  is  the 
transaction.  For  example,  each  voice  call  is  represented  in  the 
simulation  by  a separate  transaction.  The  way  a transaction  flows  in 
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the  model  is  controlled  by  a "process."  In  the  SENET  model  three 
processes  are  defined:  a control  process,  a voice  call  or  region  I 
process:  and  the  frame  process,  which  includes  the  region  II  model. 
With  this  brief  introduction  let's  look  at  the  actual  elements  of  the 
SOL -370  language  used  in  the  simulation. 

Communication  between  transactions  in  the  SOL-370  language  occurs 
through  globals.  The  types  of  SOL  globals  used  by  the  simulation  are 
INTEGER,  REAL,  and  TABLE.  INTEGER  globals  are  used  to  store  variables 
which  assume  only  whole  numeric  values  and  which  must  be  restored  for 
the  simulation  run  to  continue  after  a break  in  execution.  REAL 
variables  are  analogous  to  INTEGER  globals  in  the  restart  capability, 
except  REAL  globals  may  contain  a fractional  part  of  a numerical 
value.  These  types  of  globals  are  used  completely  at  the  user's 
di  s cretion. 

Variables  may  also  be  declared  within  a process.  These  are  common 
to  a transaction  but  may  not  be  referenced  by  any  other  transaction. 
When  the  simulation  is  restarted,  these  variables  are  restored  just  as 
are  global  variables. 

In  addition  to  SOL  global  and  local  variables,  the  simulation  uses 
TABLE'S.  SOL  tables  are  used  to  write  integer  values  on  the  log  file. 
They  do  not  store  a value.  This  historical  file  may  then  be  analyzed 
and  the  maximum  value  of  the  table,  the  time  average,  and  the  variance 
may  oe  calculated  by  a statistics  program.  Also,  two  plotting  pro- 
grams are  available  which  use  the  information  in  the  table  declaration 
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statement  and  the  log  file  to  produce  graphs  of  time  versus  table 
value.  The  value  tabulated  in  a table  Is  recorded  by  the  SOL  TABULATE 
statement. 

The  SOL  language  also  has  control  statements  that  are  used  by  the 
simulation  and  include:  WAIT,  CANCEL,  NEW  TRANSACTION  TO,  STOP,  AND 
BREAKOUT.  Each  of  these  will  be  discussed  briefly. 

At  the  beginning  of  the  simulation,  one  transaction  starts  at  the 
first  statement  of  each  process.  Additional  transactions  are  stated 
when  a NEW  TRANSACTION  TO  <label>  is  encountered.  The  symbol,  <labe1>> 
is  a statement  label  somewhere  within  the  process.  The  new  trans- 
action will  start  at  the  statement  label  while  the  present  transaction 
continues  with  the  next  statement.  All  existing  values  of  local 
variables  in  the  original  transaction  are  copied  to  the  new  trans- 
action. These  provide  initial  values  for  the  local  variables  attached 
to  the  new  transaction.  A transaction  continues  until  a transaction 
control  statement  is  executed. 

Transactions  continue  to  execute  each  statement  sequentially 
until  the  transaction  encounters  a CANCEL  statement,  the  end  of  the 
process,  or  a WAIT  statement.  If  a transaction  encounters  a CANCEL 
statement  or  the  end  of  the  process,  the  transaction  dies  and  all 
storage  used  by  its  local  variables  is  returned  to  the  system.  In 
many  cases  a transaction  must  be  initiated  by  some  event,  for  example, 
the  start  of  frame  or  the  next  data  message  arrival.  For  this  purpose 


the  SENET  simulation  and  is  one  of  the  extentions  to  the  language. 

This  is  the  BREAKOUT  TO  <label>  statement.  When  the  statement  is 
encountered  by  a transaction,  all  simulation  variables  are  output;  in 
addition,  the  label  is  output  as  the  point  where  the  present  trans- 
action will  continue  execution  when  the  simulation  is  restarted. 

After  the  BREAKOUT  statement  is  executed,  the  present  transaction 
continues  sequentially.  The  breakout  label  only  applies  to  the 
restarted  simulation  and  not  to  the  present  execution. 

In  addition  to  special  statements,  SOL  has  system  variables  and 
functions.  The  predefined  system  variables  used  in  the  simulation  are 
TIME  and  PRIORITY.  The  TIME  variable  is  the  value  of  the  simulated 
clock.  It  is  an  integer,  and  for  the  SENET  simulation  it  indicates  the 
number  of  simulated  milliseconds  which  have  elapsed  since  the  start  of 
the  simulation.  Zero  time  always  refers  to  the  start  of  the  initial 
simulation.  After  a restart,  TIME  has  the  same  value  as  when  the 
BREAKOUT  statement  was  executed.  Because  the  simulation  is  in  fact  a 
sequential  process  and  transactions  are  not  truly  run  in  parallel,  it 
is  sometime  necessary  to  ensure  that,  if  two  or  more  transactions  are 
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scheduled  to  wake  up  the  same  time,  one  will  always  be  giver  prece- 
dence over  the  others.  For  example,  the  region  I size  mu.it  oe  calcu- 
lated before  region  II.  To  ensure  that  certain  transactions  are 
executed  first,  they  may  be  assigned  precedence  by  giving  the  system 
variable  PRIORITY  a value  between  0 and  31.  When  new  transactions 
are  generated,  they  inherit  the  priority  of  the  source  transaction. 
When  a transaction  is  being  executed,  it  may  change  its  priority  at 
any  time.  A transaction  with  smaller  priority  value  will  he  executed 
before  one  with  a higher  value.  In  addition  to  system  variables,  SOL 
also  has  random  sample  generating  functions. 

The  functions  used  by  the  simulator  are  the  EXPONENTIAL! • ) 
and  PR ( - ) function.  The  EXPONENTIAL  function  returns  a real  valued 
sample  from  an  exponential  distribution  whose  mean  ^s  given  by  the 
argument.  The  PR  function  returns  a logical  true  value  if  a sample 
from  a uniform  density  on  the  unit  interval  is  less  than  or  equal  to 
its  argument. 

Because  SOL-370  is  translated  into  PL/I,  any  valid  PL/I  statement 
is  a valid  SOL  statment.  PL/ I statements  may  be  intermixed  with 
SOL-370  statements.  Blocks  of  PL/I  statements  may  be  contained 
between  special  SOL-370  statements:  PL IBEG IN  and  PLIENO.  These 
statements  have  no  effect  on  the  simulation  except  to  allow  the  con- 
tained statements  to  be  skipped  by  the  translator.  This  is  sometimes 


an  advantage  for  comment  statements  which  SOL-370  would  otherwise 
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reformat,  for  large  blocks  of  PL/I  statements  which  take  CPU  time 
for  the  translator  scan,  and  for  very  complex  PL/I  statements  which 
exceed  the  complexity  of  the  parser  in  the  SOL-370  translator.  One 
word  of  caution:  only  complete  blocks  or  DO-groups  should  be  placed 
between  the  PLI8EGIN  and  PL  I END  statements  or  the  translater  will  not 
properly  match  blocks  and  groups  in  the  overall  model. 

Variables  may  be  also  declared  using  PL/I  declare  statements; 
such  variables  are  global  to  the  entire  simulation.  If  a variable  is 
declared  by  a PL/I  statement  within  a process,  it  is  global  to  all 
transactions  executing  within  the  process.  SOL  declared  variables  are 
saved  when  a breakout  is  executed  and  restored,  and  the  simulation  is 
restarted;  PL/I  declared  variables  are  not  so  saved.  SOL  local 
variables  are  local  to  each  transaction;  PL/I  local  variables  are 
local  to  the  process.  Therefore,  PL/I  variables  declared  within  the 
process  are  common  to  all  transactions  in  the  process,  not  just 
one. 

3.  SIMULATOR  FILES 

Nine  files  are  used  to  run  the  simulator  (see  Table  II).  These 
files  provide  for  input  and  output  of  a user's  variables,  restart 
files,  and  statistics  files.  Not  all  of  these  files  are  needed  for 
any  particular  run.  Normally,  file  CARD  points  to  cards  in  the  JCL 
deck.  Files  PRINTER  and  PLIDUMP  are  set  to  dummy  because  no  useful 
information  is  written  on  them.  All  data  written  to  the  break  and  log 


files  Is  copied  to  new  files  at  restart,  permitting  magnetic  tape  to 
be  used.  After  a successful  restart  and  new  breakout,  the  old  log  and 
break  files  may  be  destroyed. 

The  JCL  needed  to  run  the  simulator  is  discussed  in  a later 
section  of  this  report.  In  the  next  subsection  the  simulator  program 
is  described  In  detail. 

4.  GLOBAL  DECLARATIONS 

This  section  discusses  the  global  declaration  portion  of  the 
simulator  (see  Figure  1).  The  simulator  uses  two  types  of  SOL  globals: 
tables  and  global  variables.  In  addition,  PL/I  globals  are  also 
used. 

The  first  statement  is  a comment  card  which  contains  within  the 
comment  text  keywords  used  by  the  SOL-370  translater.  The  keyword 
$NUM0C  Instructs  the  translater  to  number  the  translated  card  deck. 

The  SLIST  and  SSOURCE  keywords  cause  the  translator  to  produce  a 
printed  output  containing  the  translated  output  and  the  unmodified 
source  input.  These  listings  are  for  the  user  and  in  no  way  affect  the 
translator's  operation. 

The  second  statement  is  a PL/ I declare  statement.  Because  they 
are  located  in  the  global  section  of  the  program,  the  variables 
defined  within  the  statement  pertain  to  entire  program.  In  this  case, 
the  two  variables  are  defined  to  be  PL/I  intrinsic  functions.  If  this 
statement  were  removed,  there  would  be  no  functional  change  to  the 
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TABLE  II.  SIMULATOR  FILES 


- 

NAME 

REQ 

DESCRIPTION 

CARD 

A 

User  input  values 

SYSPRINT 

A 

Listings 

PRINTER 

A 

SOL  trace  table 

PLIDUMP 

A 

PL/I  error  dumps 

$LOGF 

A 

Time  history  of 
tabulated  values 

$BRKOUT 

A 

A list  of  al 1 
variables  required 
to  restart  simulation 
(produced  by  breakout 
statement ) 

SBRKFIL 

R 

A list  of  the  values 
needed  to  restart  and 
continue  the  simulation 
(copied  to  SBRKOUT 
during  restart) 

SOLDLOG 

R 

The  log  produced  by 
the  previous  run 
(copied  to  SLOGF 
during  restart) 

NOTES : 

NAMEFIL 

I 

The  list  of  table 
names  for  use  by 
statistics  and 
plot  programs 

A 

= File  is  always  required. 

R 

= File  needed 

for  restart. 

I 

= File  needed 

for  initial  run. 

13 


TALKSPURT  MODEL  WITH  SENET  RAME  8/78  SNUMOC  $LIST  SSOURCE*/ 
DCL  (MIN,  CEIL)  BUILT  IN; 

INTEGER  SIMT,  TBREAK,  TRAME,  NEXT_SOF,  VSLOT,  DLEN, 

L VOICE,  NDATA; 

REAL  VIN,  VHOLD,  DIN,  TDATA; 

TABLE  ( 0 BY  1 TO  lOOOOOO)  TQUEUE, 

( 0 BY  1 TO  1000000)  TV  RAME, 

( 0 BY  1 TO  1000000)  TSRAME, 

( 0 BY  1 TO  1000000)  TNVOICE; 

INTEGER  MQUEUE,  MS  RAME,  MNVOICE,  MV  RAME  ; 

INTEGER  QUEUE,  SRAME,  NVOICE,  VRAME  ; 


Figure  1 , Global  Declaration 
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program  except  for  a PL/I  compiler  informative  error  message.  Follow- 
ing this  statement,  the  real  “meat"  of  the  global  declaration  portion 
of  the  simulator  begins. 


The  next  statement  on  the  following  two  lines  is  the  first  SOL 
global  declaration.  The  eight  variables  are  defined  as  full  word 
integers.  Because  these  variables  are  declared  in  an  INTEGER  state- 
ment, their  values  will  be  saved  whenever  a breakpoint  is  encountered 
during  execution  of  the  simulator.  When  the  simulation  is  restarted, 
their  values  will  be  restored  before  transaction  execution  is  restarted. 

Actually,  some  of  these  variables  are  reset  by  the  user  immediately 
after  restart.  The  variables  are  described  in  more  detail  below: 

• SIMT  is  a value  input  from  file  CARD  each  simulation  start  or 
restart.  It  is  the  time  in  seconds  of  the  simulator  clock  when 
the  next  breakpoint  should  be  taken  and  the  present  simulation 
halted. 

• TBREAK  is  a value  input  from  file  CARD  each  simulation  start  or 
restart.  It  is  the  number  of  seconds  between  snapshot  print- 
outs. 

• TFRAME  is  a value  input  from  file  CARD  only  the  first  simula- 
tion start.  It  is  the  time  in  milliseconds  of  the  simulated 
SENET  frame. 

• NEXT_S0F  is  an  internally  generated  value  indicating  the  time 
of  the  next  SENET  start  of  frame  marker. 

J 
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• VSLOT  Is  a value  input  only  during  the  initial  run  of  the 


simulator.  It  represents  the  number  of  bits  in  a region  I 
slot. 

• DLEN  Is  a value  input  each  simulator  start  or  restart.  It  is 
the  mean  number  of  bits  In  a data  message. 

e LVOICE  is  an  internally  generated  counter  of  the  number  of 

r 

voice  calls  blocked  since  time  equal  zero. 

• NDATA  is  also  an  internal  variable.  It  indicates  the  number  of 
bits  allocated  for  data  packets  in  the  previous  SENET  frame, 
which  will  be  removed  from  the  data  queue  this  frame. 

The  REAL  statement  defines  its  list  of  variables  to  be  floating 
single  precision.  In  this  case,  all  the  variables  are  read-in  each 
time  the  simulation  is  initialized  or  restarted  except  the  last  one. 

The  variables  are  used  for  the  following: 

• VI N is  number  of  voice  calls  which  arrive  each  second. 

t VHOLP  is  the  mean  hold  time  in  seconds  for  a voice  call. 

t DIN  is  the  average  number  of  data  messages  per  second. 

t TDATA  Is  an  internally  generated  variable  indicating  the  time 
of  the  next  data  message  arrival. 

The  next  statement  {’ines  six  through  nine)  declares  the  tables 
used  by  the  simulator.  Tables  are  global  names  in  which  a value  may 
be  recorded.  The  first  three  tables  represent  bits  in  the  system,  and 
the  last  table,  the  number  of  voice  customers.  The  three  integers  in 
the  parenthesis  are  the  minimum,  increment,  and  maximum  values  of  the 
table.  The  values  specified  do  not  affect  the  function  of  the  table. 
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and  arbitrary  values  are  used.  The  only  required  value  is  for  the 
increment  which  must  be  nonzero.  The  following  simulated  SENET  quanti- 
ties are  placed  in  tables: 

• VFRAME  - the  number  of  bits  in  present  SENET  frame  occupied  by 
active  voice  calls 

• SFRAME  - the  number  of  bits  in  the  SENET  frame  used  by  both 
voice  and  data 

• QUEUE  - the  number  of  data  message  bits  queued 

• NVOICE  - the  number  of  voice  calls  "off-hook"  or  "in-progress." 

Tables  are  used  so  that  time  histories  of  the  values  placed  into 
the  tables  are  generated.  After  the  simulation  is  completed,  an 
independent  program  generates  the  statistics  for  these  variables. 
Because  tables  are  not  limited  to  a maximum  value,  additional  variables 
must  be  defined. 

The  next  statement  in  the  global  declaration  portion  of  the 
simulator  defines  variables  which  are  used  as  the  maximum  value  for 
the  tables.  These  maximum  values  are  read  in  during  the  initial  run  of 
the  simulator. 

The  last  line  in  this  section  of  code  defines  the  integers  used 
to  store  the  values  which  will  be  tabulated  in  the  tables  described 
above.  This  completes  the  global  declaration  portion  of  the  simulator. 
The  three  processes  are  described  in  the  following  paragraphs. 
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5.  PROCESS  SOF 

This  process  controls  the  simulation.  It  sets  the  start  of 
frame  time,  produces  snapshot  listings  of  the  table  values,  causes 
restart  variables  to  be  saved,  and  performs  all  input  to  the  simula- 
tion. The  process  is  placed  first  in  the  program  so  that  the  first 
transaction  generated  by  the  SOL  operating  system  will  start  in  this 
process.  This  transaction  will  also  be  the  first  activated  when  the 
simulation  is  restarted.  In  either  case  the  first  activities  of  the 
process  are  to  input  the  appropriate  variables  and  to  list  them.  The 
transaction  then  enters  a loop.  In  this  loop  the  transaction  calcu- 
lates the  time  of  the  next  start  of  frame,  checks  snapshot  require- 
ments, checks  for  termination  time,  and  after  termination  goes  to 
sleep  for  one  SENET  frame  period.  When  reactivated,  the  transaction 
returns  to  the  start  of  the  loop.  Because  this  process  is  always 
activated  first  and  because  it  checks  for  the  end  of  simulation,  it  is 
referred  to  as  the  control  process. 

The  first  statement  (see  Figure  2)  names  the  process,  SOF,  and 
specifies  that  only  a single  transaction  will  use  the  process  and  that 
no  resources  are  used. 

Lines  15  and  16  are  PL/I  statements  which  declare  that  variables 
N and  A are  local  to  the  process.  PL/ I statements  are  used  because 
these  variables  need  not  be  saved  and  restored  between  simulation 
runs.  The  next  statement  sets  the  priority  of  the  process.  The 
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12  PROCESS  SOF  T=  1 R=0' 

13  /*  CONTROL  PROCESS  ’ INPUTS  DATA,  STOPS  SIMULATION, 

14  PRINTS  INTERMEDIATE  RESULTS,  SETS  SOF  TIME  */ 

15  DCL  N FIXED  8 INARY (15)  STATIC  INIT(O); 

16  DCL  A FLOAT(21)  BINARY  STATIC  INIT(0.); 

17  PRIORITY  = 1;/*  DO  THIS  PROCESS  FIRST  */ 

18  READIN: 

19  PLIBEGIN; 

20  DCL  FLOAT  BUILT  IN; 

21  GET  FILE  (CARD)  LIST(A)  ; 

22  SIMT  = CEIL(1000.*A);  /*  CHANGE  TO  MS  */ 

23  GET  FILE  (CARD)  LIST(A)  ; 

24  T8REAK  = CEIL(100Q.*  A );  /*  CHANGE  TO  MS  */ 

25  GET  FILE  (CARD)  LIST(A)  SKIP; 

26  V IN  = 1000.  / A;  /*  CHANGE  TO  MS  V 

27  GET  FILE  (CARD)  LIST(A)  ; 

28  YHOLD  = 1000.  * A;  /*  CHANGE  TO  MS  */ 

29  GET  FILE  (CARD)  LIST(A)  SKIP; 

30  DIN  = 1000.  / A;  /*  CHANGE  TO  MS  */ 

31  GET  FILE  (CARD)  L 1ST ( DLEN)  ; 


32 

IF  TIME  * 0 THEN 

33 

GET  FILE(CARD)  SKIP(l) 

34 

L 1ST ( MQUEUE , MS  FRAME ,MV  FRAME ,MNV0 ICE , VSLOT ,T  FRAME ) ; 

35 

PUT  FILE(SYSPRINT) 

EDIT ('SIMULATION  TIME 
'PRINT  INTERVAL 

36 

' ,.001*FL0AT(SIMT) , ' 

SEC', 

37 

' ,.001*F10AT(TBREAK),' 

SEC', 

38 

'VOICE  ARRIVAL  RATE 

' , 1000. /V  IN  , ' 

PER  SEC' 

39 

'DATA  ARRIVAL  RATE 

',1000. /DIN  , ' 

PER  SEC' 

40 

'VOICE  HOLDING  TIME 

' ,.001*VH0LD  , ' 

SEC', 

41 

'DATA  MESSAGE  LENGTH 

' ,DLEN, 

BITS’ , 

42 

'VOICE  SLOT  SIZE 

'.VSLOT, 

BITS', 

43 

'FRAME  TIME 

' , .001*FL0AT( TFRAME) , ' 

SEC' , 

44 

'MAX  QUEUE  SIZE 

' .MQUEUE 

8 ITS ' , 

45 

'MAX  VOICE  FRAME 

' , MV FRAME 

BITS' , 

46 

'MAX  TOTAL  FRAME 

' , MS  FRAME 

BITS', 

47 

'MAX  NUMBER  OF  VOICE 

' .MNVOICE 

CALLS ' ) 

48 

(SKIP,(2)(A(22) ,F(12, 

3) ,A(14) ) ) ; 

’ LVO ICE ' , ' QUE', 

V FRAME ‘ , ' TIME') 

49 

PUT  FILE(SYSPRINT)  SKIP  ED  IT ( 
' NVOICE','  S FRAME ' , ' 

50 

51 

(X(13) , (6)A(10) ); 

52 

PLIEND; 

53 

START: 

54 

NEXT  SOF  = TIME  + TFRAME; 

55 

IF  N >=  TBREAK  THEN 

56 

DO; 

57 

PUT  FILE(SYSPRINT)  SKIP 

58 

ED  IT ( L VO  ICE .QUEUE ,NVO ICE , S FRAME , V FRAME ,T IME-T FRAME ) 

59 

( X ( 10 ) ,(6)F(10) ); 

60 

N = 0; 

61 

END; 

62 

N = N + TFRAME; 

63 

IF  TIME  > SIMT  THEN  /*  END  THIS  RUN  */ 

64 

DO; 

65 

BREAKOUT  TO  READ  IN; 

66 

STOP; 

67 

END; 

68 

WAIT  TFRAME; 

69 

GO  TO  START; 

70 

END; 

71 

Figure  2 (contd).  Process  SOF 
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priorities  of  the  other  two  processes  have  higher  values.  Therefore, 
if  transaction  reactivations  are  scheduled  for  identical  times,  the 
transaction  in  this  process  will  be  reactivated  first. 

The  label  on  the  next  line  is  used  to  indicate  the  restart  point 
for  this  transaction  when  the  simulation  is  restarted.  This  is 
explained  in  more  detail  later. 

The  next  line  is  an  instruction  to  the  translator  that  all  text 
between  PLIBEGIN  at  line  19  and  PLIEND  at  line  52  may  be  skipped  by 
the  translator  and  immediately  copied  to  the  translated  output.  This 
action  speeds  up  translation  and  stops  the  translator  from  reformat- 
ting the  text.  All  the  statements  contained  between  lines  19  and  52 
are  complete,  valid  PL/I  statements.  No  SOL  statements  are  mixed 
in. 

These  statements  between  lines  19  and  52  perform  all  user  inputs 
and  outputs  of  the  simulator.  Line  20  is  a declaration  of  a PL/I 
built-in  function.  Statements  on  lines  21  through  31  input  the 
variables,  read  each  simulation  start  or  restart,  and  convert  the 
values  to  internal  units.  Lines  32  through  34  input  values  which  are 
read  only  once  at  the  beginning  of  the  simulation.  These  are  the 
maximum  values  of  the  SENET  frame  structure.  Lines  35  through  51  list 
all  input  variables  and  write  the  heading  for  the  snapshots. 

The  first  set  of  variables  are  input  in  units  of  seconds  or 
arrivals/  second.  Variables  input  in  seconds  are  converted  to  milli- 
seconds, and  arrivals/second  are  converted  to  milliseconds  between 
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arrivals.  The  SKIP  option  on  the  GET  statements  causes  a skip  to  the 
next  card;  that  is,  two  values  are  input  per  card.  Text  or  card 
numbers  may  follow  the  two  inputs.  The  last  input  card  is  an  exception. 

The  fourth  input  card  is  only  input  at  simulator  time  zero.  Six 
integer  values  are  input:  the  first  four  establish  table  maximums, 
and  the  last  two  set  the  voice  slot  size  in  bits  and  the  SENET  frame 
duration  in  milliseconds.  No  error  checking  is  done  for  these  vari- 
ables, and  the  uSer  must  be  careful  in  assigning  values. 

All  values  input  on  the  fourth  card  must  be  positive  integers. 
MVFRAME  represents  the  number  of  bits  of  the  SENET  frame  which  may  be 
allocated  for  region  I data.  MSFRAME  represents  the  total  number  of 
bits  in  a frame.  Therefore,  MVFRAME  may  not  be  greater  than  MSFRAME 
(i.e. , the  simulator  cannot  handle  TASI  transmission).  Thus,  MNVOICE, 
the  maximum  number  of  voice  channels,  may  not  exceed  the  region  I 
capacity.  Therefore, 

MNVOICE*VSLOT<=MVFRAME, 

must  be  satisfied. 

After  the  variables  are  input,  the  transaction  enters  the  loop  in 
lines  53  to  70.  This  loop  will  be  executed  until  the  logical  test 
at  line  64  becomes  true.  When  the  logical  test  become  true,  the 
breakout  statement  causes  all  variables  and  the  transaction  status  to 
be  saved  in  the  break  file.  In  addition,  label  READIN  is  established 
as  the  restart  point.  When  the  simulation  is  restarted,  the  break 
file  is  read  and  simulation  continues  with  this  transaction  at  the 
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statement  pointed  to  by  label  READ  IN.  The  breakout  statement  does  not 
halt  the  present  simulation.  The  transaction  continues  execution  to 
line  67  where  the  STOP  statement  terminates  the  simulation. 

The  other  statements  are  standard  PL/I.  The  first  line  sets  the 
next  start  of  frame  time  global  which  is  used  to  synchronize  all  the 
transactions.  The  next  10  lines  form  a primitive  modulo  counter  for 
the  snapshot  listing. 

6.  PROCESS  VOICE 

This  process  handles  the  region  I allocation  of  the  SENET  frame. 
This  is  combined  with  the  region  II  traffic  in  process  FRAME.  Each 
call  is  handled  as  an  independent  transaction  within  this  process.  In 
theory,  as  each  call  arrives  a new  transaction  is  initiated  to  await 
the  next  call  arrival,  and  this  call  is  added  to  the  number  of  "in 
prog-ess"  calls.  The  present  transaction  executes  the  speech  activity 
model  (see  Appendix  A)  for  this  call.  This  is  continued  until  it  is 
time  for  the  call  to  "hang  up".  At  this  time  the  transaction  removes 
the  call  from  the  "in  progress"  calls  and  cancels  itself.  The  actual 
simulation  process  varies  slightly  from  this  in  order  to  decrease  CPU 
time. 

This  process  is  the  most  complex  in  the  simulation.  Its  discus- 
sion is  therefore  divided  into  three  parts:  preamble,  speech  simula- 
tor, and  conclusion.  In  the  discussion  of  the  preamble  the  meaning  of 
the  local  variables,  call  arrival  and  allocation,  are  detailed.  The 


72  PROCESS  VOICE,  T=195,R«0; 

73  REAL  TM; 

74  INTEGER  TI.STATE, TALK, LTALK.ACTIYE, TEND; 

75  DCL  T aOAT(21)  BINARY  INIT(O)  STATIC; 

76  PRIORITY  - 2 ;/*  DO  THIS  PROCESS  AFTER  FRAME  START*/ 

77  START: 

78  WAIT  EXPONENTIAL  (VIN);/*  WAIT  FOR  CALL  ARRIVAL  */ 

79  NEW  TRANSACTION  TO  START;/*  WAIT  FOR  NEXT  ARRIVAL*/ 

80  /*  PROCESS  THIS  CALL  */ 

81  IF  NVOICE  >=■  MNVOICE  THEN  /*  NO  MORE  SPACE  IN  FRAME*/ 

82  00* 

83  ’ LVO ICE  * L VO  ICE  + 1; 

84  GO  TO  FINIS; 

85  END; 

86  /*AOD  THIS  CALL  TO  TOTAL  CALLS  IN  PROGRESS*/ 

87  NVOICE  « 1 + NVOICE  ; 

88  TABULATE  NVOICE  IN  TNVOICE; 

89  /*  INITIALIZE  FOR  THIS  CALL  */ 

90  TALK  - 0;  /*  NO  TALKSPURT  IS  IN  SENET  FRAME  */ 

91  STATE,  ACTIVE  * 1;  /*SET  INITIAL  STATE  */ 

92  /*  STORE  CALL  START  TIME  FLOATED  INTO  VARIABLE  TI  */ 

93  TM  - TIME; 

94  TEND  * EXPONENT IAL ( V HOL D ) + TIME;  /*  CALL  END  TIME  */ 

95  /*  WAIT  TIL  START  OF  NEXT  FRAME  TO  RUN  MODEL  */ 

96  WAIT  NEXT  SOF  - TIME; 


speech  simulator  is  then  briefly  discussed.  Finally,  the  region  I 
allocation,  the  transaction  "sleep"  period,  and  call  hang  up  are 
di scussed. 

The  discussion  will  start  with  the  preamble  (see  figure  3). 
Because  each  call  is  handled  by  a separate  transaction,  the  number 
assigned  to  T in  the  PROCESS  statement  must  satisfy  the  relation: 

KMNVOICE+l. 

This  allocates  a sufficient  number  of  transactions  to  this  process 
to  simulate  MNVOICE  simultaneous  voice  calls,  a transaction  which 
awaitsthe  next  call,  and  the  present  transaction.  The  number  of 
simultaneous  calls  which  may  be  simulated  by  the  present  simulator  is 


Lines  73  through  75  declare  the  local  variables  to  this  process. 
The  variables  declared  in  the  REAL  and  INTEGER  statements  are  local  to 
each  transaction  and  are  used  to  store  the  "state"  of  the  speech 
simulator  for  this  transaction: 


• TM  - The  time  of  the  next  state  transition  in  the  speech 
activi ty  model . 

• STATE  - The  next  state  to  be  transitioned  to  in  the  speech 
model . 

• TALK  - Non  zero  if  space  is  to  be  allocated  in  region  I of 
this  frame  for  this  call. 

• LTALK  - Non  zero  if  space  were  allocated  in  region  I of  the 
last  frame  for  this  call. 


30 


• ACTIVE  - Non  zero  indicates  that  the  present  state  of  the 
speech  model  represents  a talkspurt.* 

e TEND  - The  time  the  call  will  hang  up. 

* 

In  addition  to  the  above  integer  variables  the  process  uses  a tempor- 
ary floating  point  variable: 

e T - The  time  in  floating  seconds  until  the  next  state  transi- 
tion in  the  speech  activity  model. 

The  last  declare  statement  is  again  a PL/I  built-in  function  specifica- 
tion to  suppress  a compiler  error  message. 

Line  76  is  only  executed  at  time  zero  and  sets  the  priority  of  all 
transactions  in  this  process.  The  value  assigned  ensures  that  region 
I allocations  are  complete  before  the  region  II  allocations  are 
calculated  and  the  SENET  frame  is  constructed. 

In  this  model  voice  calls  are  assumed  to  have  Poisson  arrivals. 
Therefore,  they  have  an  exponential  interarrival  distribution  of  mean 
VIN.  This  is  simulated  by  lines  77  through  79. 

After  a call  arrival,  the  transaction  checks  for  space  availabil- 
ity. If  space  is  not  available  the  number  of  lost  calls  is  incremented 
and  the  transaction  cancels  itself.  This  is  accomplished  by  the 
code  starting  at  line  81. 

When  the  transaction  reaches  line  86,  it  is  known  that  there  is 
room  for  the  call.  The  call  is  added  to  NVOICE,  which  indicates  the 
number  of  calls  in  progress  (line  87),  and  the  speech  activity  simula- 
tor is  initialized  by  lines  90  through  96.  The  present  time  is  converted 
Defined  in  Appendix  A. 
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to  a floating  value  and  stored  in  TM  at  line  93.  Line  94  calcu- 
lates the  time  for  the  call  to  end.  Notice  that  the  call  holding  time 
is  simulated  to  have  an  exponential  distribution  of  mean  VHOLD.  The 
final  statement  causes  the  transaction  to  pause  until  the  next  start 

I of  frame  before  starting  the  voice  activity  portion  of  the  process. 

When  a transaction  reaches  this  point  (line  96)  in  the  process, 
the  initialization  or  preamble  is  finished.  From  here,  the  trans- 
action starts  the  voice  activity  and  region  I allocation  portion  of 
the  process. 

The  speech  activity  model  is  shown  in  Figure  4.  The  boxes  repre- 
sent  the  "state"  of  the  speaker.  The  speaker  is  uttering  sounds  in  / 

the  talking  state,  state  1.  The  speaker  then  transitions  to  one  of 
the  two  silent  states.*  After  staying  in  one  of  these  silent  states 
for  a period,  the  speaker  returns  to  state  1 and  the  process  repeats 
i tsel f . 

In  the  simulation,  the  period  of  time  a speaker  remains  in  any 
state  is  a sample  from  a random  variable  with  an  exponential  distribu- 
tion. The  mean  of  the  distribution  is  indicated  in  a corner  of  the 
box  representing  each  state.  The  numbers  near  the  transition  lines  in 
the  figure  indicate  the  rate  at  which  the  given  transition  will  be 
taken. 

In  the  simulator,  the  speech  activity  model  starts  at  line  97  (see 
Figure  5).  The  previous  state  information  is  saved  in  lines  101  and 
102.  The  model  is  run  until  the  the  next  transition  occurs  after  the 
present  start  of  frame  (the  value  of  TIME). 

* Defined  in  Appendix  A. 
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LOOP: 

IF  TIME  > TEND  THEN/*CALL  ENDED*/ 

GO  TO  HANG  UP; 

PLIBEGIN; 

LTALK  = TALK ; /*REMEMBER  IF  ACTIVE  LAST  FRAME  */ 

TALK  * ACTIVE ;/*REMEMBER  ACTIVITY  OF  PRESENT  STATE*/ 

/*  RUN  MODEL  UP  TO  PRESENT  TIME.  TI  (ALIAS  TM ) INDICATES  TIME 
OF  NEXT  MODEL  STATE  CHANGE  */ 

DO  WHILE(TIME  > TM); 

SELECT ( STATE ) ; /*  CASE  STATEMENT*/ 

WHEN  ( 1 ) 

DO;  /*  STATE  1 HAS  TWO  EXITS  */ 

T = EXPONENTIAL  (0.05); 

/*  TALK  GETS  SET  TO  1 IF  STATE  1 IS  EVER  PASSED  THRU  RUNNING  MODEL 

UP  TO  TIME  “TIME"*/ 

TALK, ACTIVE  = 1; 

PLIEND; 

IF  PR ( 0 . 03 1 5 ) THEN  /*  EXIT  TO  STATE  2 */ 
STATE  =2; 

ELSE  /*  EXIT  TO  STATE  3 */ 

STA1E  = 3 ; 

PLIBEGIN; 

END;  /*  STATE  = 1 */ 

WHEN  ( 2 ) /*  EXIT  TO  STATE  1*/ 

DO; 

ACTIVE  = 0; 

STATE  = 1; 

T = EXPONENTIAL  (2.8457); 

END;  /*  STATE  = 2*/ 

WHEN  ( 3 ) /*  EXIT  TO  STATE  1 */ 

DO; 

ACTIVE  = U; 

T = EXPONENTIAL  ( .01  ) ; 

STATE  = 1; 

END;  /*  STATE  = 3*/ 

OTHERWISE;/*  DO  NOTHING  */ 

END;  /*  SELFCT*/ 

TM  = TM  + 1000.  * T; 

END ; /*T IME  > TM*/ 


Figure  5.  Process  Voice  : Simulator 
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The  SENET  statement  in  conjunction  with  the  WHEN  statements 
determine  which  block  of  code  to  execute,  one  block  for  each  state. 

If  the  model  is  in  state  1 before  the  model  starts  to  run,  passes 
through  it  while  running,  or  ends  in  It,  then  a voice  slot  must  be 
coded  in  the  SENET  frame.  This  is  indicated  by  variable  talk  being 
set  to  "l".  This  will  be  checked  in  the  next  part  of  the  process. 

The  function  of  variables  T and  TM  have  been  explained  previously. 

One  additional  statement  should  be  explained:  line  114. 

The  talkspurt  state*  has  two  exits.  The  time  interval  that  the 
process  remains  in  the  state  has  an  exponential  density.  After  this 
interval  the  talking  state  will  transition  to  one  of  the  other  states. 
The  state  transitioned  to  depends  on  the  rate  of  transition.  Given 
that  the  state  is  left  20  times  a second,  it  transitions  to  state  3 
with  probability  19.37/20  and  to  state  2 otherwise.  The  notation  PR 
(0.0315)  indicates  a sample  drawn  from  a uniform  density  defined  on 
the  unit  interval.  When  this  sample  is  less  than  0.0315,  the  next 
state  is  chosen  to  be  2;  otherwise  3. 

When  variable  TM  becomes  greater  than  the  present  time,  the  region 
I allocation  is  executed  by  the  transaction.  (The  simulation  time  is 
at  a SENET  start  of  frame.)  The  talkspurt  model  represents  the 
activity  events  of  this  simulated  call  during  the  last  SENET  frame. 

The  simulot  is  now  ready  for  the  region  I of  the  SENET  frame  to  be 
allocated. 

* Defined  in  Appendix  A„ 
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In  the  conclusion  of  process  VOICE  (see  Figure  6),  this  call  is 
allocated  in  region  I.  In  addition,  this  section  contains  the  call 
termination  statements  plus  code  to  decrease  simulation  time. 

Variable  TALK  is  nonzero  if  this  call  needs  space  allocated  in 
this  frame.  Variable  LTALK  indicates  whether  space  was  allocated  in 
the  previous  frame.  The  allocation  for  the  call  need  only  be  changed 
if  it  has  changed  from  the  previous  frame.  This  is  done  by  lines  138 
through  145.  After  this  action,  the  transaction  goes  to  sleep  until 
the  next  frame.  At  this  point,  we  could  loop  back. 

If  we  loop  back  to  label  LOOP  after  the  wait  at  line  146,  the 
simulation  would  work  properly,  but  it  would  use  an  excessive  amount 
of  CPU  time.  This  is  because  it  takes  a lot  of  time  to  put  a trans- 
action to  sleep  and  reawaken  it.  If  we  consider  the  frame  size  of  a 
typical  SENET  simulation  and  the  speech  model,  a way  to  decrease  CPU 
time  becomes  apparent.  Typically,  the  SENET  frame  has  a length  in  the 
order  of  10  ms.  The  average  talkspurt  is  50  ms,  and  the  long  silence 
is  almost  3 seconds.  These  two  states  are  occupied  over  90%  of  the 
time.  Therefore,  the  call  occupancy  holds  constant  for  many  SENET 
frame  periods.  The  transaction  waits  a single  frame.  If  there  is  no 
state  transition  scheduled  for  this  frame,  the  activity  of  the  present 
state  is  used  to  determine  the  occupancy  of  the  SENET  frame  until  the 
start  of  frame  following  the  next  transition  time. 
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/♦CHANGE  FRAME  ONLY  IF  SPEECH  ACTIVITY  HAS  CHANGED  */ 

PLIEND; 

IF  LTALK*1 » TALK  THEN 
D0‘ 

IF  TALK  * 0 THEN  /*  HAVE  GONE  S ILENT*/ 

V FRAME  « VFRAME-VSLOT ; 

ELSE  /*  HAVE  STARTED  TALKING*/ 

V FRAME  ■ VFRAME+VSLOT ; 

TABULATE  V FRAME  IN  TVFRAME; 

END; 

WAIT  T FRAME; 

/♦CODE  ROM  HERE  ON  ADDED  TO  DECREASE  SIMULATION  TIME  */ 

IF  TIME  > TM  THEN  /*  NO  SAVINGS  POSSIBLE  */ 

GO  TO  LOOP; 

IF  TIME  > TENO  THEN 
GO  TO  HANG  UP; 

IF  TALK-* - ACTIVE  THEN 
DO; 

IF  ACTIVE  * 1 THEN 

VRAME  * VRAME+VSLOT; 

ELSE 

VRAME  * VRAME-VSLOT; 

TABULATE  VRAME  IN  TVRAME; 

ENO; 

TALK  * ACTIVE; 

IF  TEND  < TM  THEN  /*  CALL  ENOS  BEFORE  NEXT  STATE  CHANGE 
DO; 

WAIT  TEND  - TIME;/*  SLEEP  UNTIL  CALL  ENOS  */ 

GO  TO  HANG  UP; 


ELSE 


END; 

DO; 


HANG  UP: 


END; 


WAIT  TM-TIME; 
WAIT  NEXT_SOF 
GO  TO  100F; 


TIME; 


FINIS: 

END; 


WAIT  NEXT_SOF  - TIME;/*  HANGUP  IN  NEXT  RAME*/ 
NVOICE  * NVO ICE-1; 

TABULATE  NVOICE  IN  TNVOICE; 

IF  TALK  = 1 THEN  /*  TALKSPURT  IN  PROGRESS*/ 

DO’ 

* VRAME  = VRAME-VSLOT; 

TABULATE  VRAME  IN  TVRAME; 

END; 


Figure  6.  Process  Voice  : Conclusion 
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The  final  portion  of  the  process  terminates  the  call.  The  present 
frame  allocation  is  not  changed.  The  transaction  waits  until  the 
next  frame  and  then  removes  the  call  from  the  number  of  calls  in 
progress  and  removes  any  region  I allocation. 

All  transactions  in  this  process  update  the  region  I allocation, 
VFRAME,  before  the  final  SENET  frame  occupancy  is  determined  in 
process  FRAME. 

7.  PROCESS  FRAME 

The  last  process  (see  Figure  7)  accumulates  incoming  messages  in 
the  queue  and  adds  the  region  II  traffic  from  this  queue  to  the  region 
I traffic  allocated  in  process  VOICE  forming  the  final  SENET  frame 
allocation.  Both  of  these  functions  are  performed  by  a single  trans- 
action. The  transaction  loops  in  lines  189  through  197  accumulating 
message  bits  in  the  queue  for  all  messages  which  arrive  before  the 
next  start  of  frame  and  updating  the  time  of  next  message,  TDATA. 

When  the  next  message  is  scheduled  to  arrive  after  the  start  of  frame, 
the  transaction  drops  out  of  the  loop  and  waits  for  the  start  of 
frame.  This  permits  the  transactions  in  the  voice  process  to  complete 
the  region  I allocation  and  for  a new  start  of  frame  time  to  be 
calcul ated. 

At  this  time  the  transaction  removes  the  bits  transmitted  in 
region  II  during  the  previous  frame,  NDATA,  from  the  queue.  A new 
NDATA  is  calculated  for  this  frame  based  on  the  number  of  bits  remain- 
ing in  the  SENET  frame  for  region  II  traffic  and  the  number  queued. 


38 


184  PROCESS  FRAME ,T*1,R»0; 

185  DCL 

186  DELTA  FIXED  BINARY (31)  STATIC  INIT(O); 

187  PRIORITY  =3;/*  ACTIVATED  AFTER  CLASS  I FINISHED  */ 

188  TDAT A*EXPONENT IAL  { D IN ) ; /*  TIME  OF  FIRST  MESSAGE*/ 

189  START: 

190  IF  TDATA  < NEXT  SOF  THEN 

191  DO; 

192  WAIT  TDATA  - TIME; 

193  QUEUE  = QUEUE+EXPON£NTIAL(DLEN);/*QUEUE  IT*/ 

194  TABULATE  QUEUE  IN  TQUEUE; 

195  TDATA  » EXPONENT IAL (DIN)  + TDATA; /*TOATA  OF  NEXT  MSG*/ 

196  GO  TO  START;/*  WAIT  FOR  NEXT  ONE*/ 

197  END; 

198  /*THIS  SECTION  DETERMINES  THE  SENET  FRAME  */ 

199  WAIT  NEXT_SOF-TIME;  /*WAIT  FOR  NEXT  FRAME  */ 

200  IF  NDATA  > 0 THEN 

201  DO;  /*REM0VE  DATA  PUT  IN  PREVIOUS  FRAME  */ 

202  QUEUE  = QUEUE-NDATA; 

203  TABULATE  QUEUE  IN  TQUEUE; 

204  END; 

205  PLIBEGIN; 

206  ^CALCULATE  SPACE  NEEDED  OR  AVAILABLE  THIS  FRAME*/ 

207  NDATA=MIN(QUEUE,MSFRAME-VFRAME); 

208  /*CHANGE  FROM  LAST  FRAME  */ 

209  /*  VFRAME  VOICE  THIS  FRAME  */ 

210  /*  S FRAME  LAST  FRAME  OCCUPANCY*/ 

211  /*  NDATA  THIS  FRAME  DATA*/ 

212  DELTA  = NDATA  + VFRAME  - S FRAME; 

213  PLIEND; 

214  IF  DELTAS  = 0 THEN 

215  DO; 

216  S FRAME  = SFRAME+DELTA;  /*CHANGE  SENET  FRAME*/ 

217  TABULATE  S FRAME  IN  TS FRAME; 

218  END; 

219  GO  TO  START; 

220  END; 


Figure  7.  Process  Frame 
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The  change  in  allocation  between  this  SENET  frame  and  the  previous 
frame  is  calculated  and  stored  in  temporary  variable,  DELTA.  This 
value,  if  not  zero,  is  used  to  update  the  frame. 

Because  of  the  wait  until  the  present  start  of  frame  and  the  low 
priority  assigned  to  this  transaction,  the  transaction  in  the  first 
process  updates  the  start  of  frame.  Therefore,  when  this  transaction 
loops  back  to  label  START,  TDATA  is  compared  to  the  next  start  of 
frame  and  the  process  continues. 


IV.  RUNNING  THE  SIMULATOR 


Examples  of  running  the  simulator  are  presented  in  the  first  part 
of  this  section.  These  are  followed  by  examples  for  running  the 
statistics  and  plot  programs.  The  section  concludes  with  an  example 
of  translating  and  compiling  the  model.  After  the  extended  translator 
is  formally  released,  some  of  the  dataset  names  may  need  to  be  changed. 

Due  to  the  large  amount  of  data  recorded  on  the  log  file,  runs 
for  long  simulation  periods  should  use  magnetic  tape.  Many  times  I 
failed  to  give  enough  CPU  time  for  a simulation.  Therefore,  the  job 
ended  without  properly  closing  the  files.  If  the  log  file  were  a 
member  of  a disk  dataset,  all  results  of  the  simulation  would  be 
lost.  Using  magnetic  tape  stops  this  loss.  Even  if  the  simulation 
aborts,  the  statistics  or  plot  programs  are  able  to  recover  most  of 
the  data. 

However,  magnetic  tapes  can  be  unwieldy  both  from  the  user's  and 
the  computer  operator's  viewpoints.  Therefore,  a compromise  is  used. 
The  total  simulation  usually  takes  two  runs:  an  initialization  run  to 
load  the  SENET  frame  to  near  steady  state,  and  then  a run  using  the 
variables  of  interest.  The  author's  compromise  is  to  run  the  initial 
simulation  using  all  disk  files.  The  following  run  then  uses  these 
disk  files  as  input  and  produces  new  files  on  tape.  The  tapes  are 
used  as  Input  to  the  statistics  program. 

The  JCL  decks  for  this  type  of  operation  are  shown  in  Figures  8 
and  10;  output  listings  are  shown  in  Figures  9 and  11.  Notice  that 
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//R3363S0L  JOB  ( A1 00 , , 50 , 9 , 220  ) , ' C ASTELLO ' , 

//  NQTIFY=R3363,MSGCLASS=Q 

//GO  EXEC  PGM=M0DEL,REGI0N=220K,TIME=(0,50) ,PARM='/0' , 

//  CQNO= ( EVEN ) 

//STEPLIB  DO  DSN=R3363.L0AD,DISP=SHR 
//SYSPRINT  DO  SYSOUT  = ( A ,U  ) 

//PRINTER  DD  SYSOUT=(A,U) 

//NAMEFIL  DD  DS N=R3363 . SOL . NAMES ( MODEL ), D ISP =S HR 
//PLIDUMP  DD  DUMMY 

//SLOGF  DD  DSN  = R3363.S0L.L0G(M0DEL0)  , D ISP  = S HR  OUT 

//JBRKOUT  DD  DSN=R3363. SOL . RESTART (MODELO ), DISP=SHR 
//CARD  DD  * 

1.000,  .990  END  SIMULATION  TIME, PRINT  INTERVAL 

150.000. 300,  VOICE  ARRIVAL  RATE,  CALL  HOLDING  TIME 

10.000. 1500,  DATA  ARRIVAL  RATE,  LENGTH 
1000000,15440,15440,193,  80,  10 

r 


Figure  8.  JCL  For  Simulator  Run  Using  Disk 
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//R3363S0L  JOB  ( AlOO , ,450 , 9 , 230 ) , ' CASTELLO ' , 

//  N0TIFY=R3363,MSGCLASS=Q 

/♦SETUP  NEED  TWO  TAPES  RING  IN  2893  2814 

/ /G02  EXEC  PGM=M0DEL,TIME=(6,50) , 

//  PARM='/1' ,REGI0N=230K 

//STEPLIB  DD  DSN=R3363.L0AD,DISP=SHR 

//SYSPRINT  DD  S YSOUT= ( A , U ) 

//PRINTER  DD  SYSOUT=(A,U) 

//SOLDLOG  DD  DSN=R3363 . SOL . LOG ( MODELO ) , D ISP=S HR 
//SBRKFIL  DD  DSN  = R3363 . SOL . RESTART ( MODELO ), D ISP  = S HR 
//SBRKOUT  DD  D ISP= ( NEW, KEEP ), DSN  = UR3363 . REG02  , 

//  UNIT=TAPE9 ,DCB=(RECFM=FB,LRECL=133,BLKS I ZE= 1 330 ) , 
//  VOL=SER=002814,LABEL=1 
//CARD  DD  * 

30.000,  1.,  END  SIMULATION  TIME, PRINT  INTERVAL 
0.520,300,  VOICE  ARRIVAL  RATE,  CALL  HOLDING  TIME 

500.000. 1500,  DATA  ARRIVAL  RATE,  LENGTH 
/* 

//PLIDUMP  DD  DUMMY 

//JLOGF  DD  DISP=(NEW,KEEP),DSN  = l!R3363.L0G2, 

//  UNIT=TAPE62,LABEL=1,V0L=SER=002893 

//STAT2  EXEC  PGM=S0LSTAT,REGI0N=120K 
//STEPLIB  DD  DSN=R3363.L0AD,DISP=SHR 
//SYSPRINT  DD  SYSOUT=(A,N) 

//NAMME  DD  DSN  = R3363 . SOL . NAMES ( MODEL ), D ISP  = S HR 
//LOG  DD  DSN=UR3363.L0G2,DISP=SHR,UNIT=TAPE62, 

//  LABEL=1 ,V0L=SER=002893 
//PLIDUMP  DD  DUMMY 
//SYS  IN  DD  * 

5000  50000  10000 
TQUEUE 
TV  FRAME 
TS  FRAME 
TNVOICE 
/♦ 


Figure  10.  JCL  For  Simulator  Run  : Disk  In,  Tape  Out 
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the  PARM  option  on  the  EXEC  card  is  set  to  zero  for  the  initial 
simulation  and  to  one  for  the  continuation.  In  general , the  value  is 
set  to  the  number  of  the  breakout  on  the  break  file  at  which  the 
simulation  is  to  be  restarted.  If  the  simulation  were  to  be  continued 
after  the  second  run,  this  field  would  be  set  to  three. 

The  time  parameter  is  required  on  the  EXEC  card  for  execution  time 
over  5 minutes.  The  time  required  to  execute  the  example  program  was 
23  seconds  for  the  initial  run  and  390  seconds  for  the  continuation. 

To  estimate  run  times  for  other  cases,  assume  that  run  time  is  propor- 
tional to  the  voice  Erlang  loading  and  the  data  message  arrival 
rate. 

The  statistics  program  uses  the  NAMME  file  created  by  the  initial 
simulation  and  the  $L0GF  from  any  simulation  run.  The  maximum, 
average,  and  variance  may  be  calculated  for  up  to  four  logged  vari- 
ables for  any  time  period  on  the  log.  The  statistics  are  printed  at  a 
user  specified  interval.  This  running  average  allows  the  user  to 
ascertain  that  the  simulation  has  been  run  long  enough  for  the  statis- 
tics to  be  nearing  steady  state. 

The  input  variables  are  shown  in  Figure  12.  The  statistics  will 
be  calculated  from  "time  begin"  to  "time  end."  Statistics  are  printed 
each  "print  increment"  and  at  the  end  of  file  on  the  log.  The  JCL 
for  a statistics  run  is  shown  in  Figure  13. 

A plot  program  is  also  available  which  produces  a time  plot  of  up 
to  four  logged  variables  on  the  CALCOMP  flatbed  plotter.  Its  JCL 
and  plot  are  shown  in  Figures  14  and  15.  The  input  variables  are 
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Card  1 -n* 

Time  begin  nm 

Card  2s- 5 

Statistics  Program  Input 


Time  end 


Print  Increment 


Figure  12.  Statistics  Program  Input 


//R3363STA  JOB  (A100, ,IZ0, 9,120 ), 'CAS TELLO  , 

//  NOT IFY*R3363 ,MSGCLASS®Q 

//STAT1  EXEC  PGM®SOLSTAT ,REG I0N*120K 
//STEPLIB  DO  QSN=R3363 . LOAD , D ISP®S  HR 
//SYSPRINT  DO  SYS0L7®(A,N) 

//NAMME  DO  QSN*R3363. SOL. NAMES (MODEL) ,DISP=SHR 
//LOG  00  DSN®UR3363 . L0G1 , 0 ISP=S HR , UN IT=TAPE62 , 
//  LABEL=1 ,VQL*SER®003707 
//PLIDUMP  DO  DUMMY 
//SYS  IN  DO  * 

5000  50000  1000 
TQUEUE 
TV FRAME 
TS FRAME 
TMVOICE 


Figure  13.  JCL  For  Statistics  Program 
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* CASTELLO  * , 


//R3363STA  JOB  ( AlOO , , 120 ,9 , 180 ) , 

//  NOT IFY=R3363 ,MSGCLASS=Q 
//STAT2  EXEC  PGM*SOLPLOT ,REG  I0N  = 180K 
//STEPLIB  DO  DSN=R3363.L0AD,DISP=SHR 
//SYSPRINT  DD  SYSOUT*(A,M) 

//NAMME  DD  DSN=R3363 .SOL . NAMES (MODEL ), D ISP  = S HR 
//LOG  DO  DSN=UR3363.L0G2,DI$P=SHR,UNIT=TAPE62, 

//  LABEL=1,V0L=SER=002893 

//PLIDUMP  DD  DUMMY 
//FT06F001  DD  SYSOUT=(A,U) 

//FT10F001  DD  0SN  = U1JCA.  P3363.  PLOT,  DISP=»(MEW,  KEEP)  , 

//  UNIT=TAPE9,LABEL=(1,NL),V0L*SER*001549, 

//  DCB= ( REC  FM=VS ,LRECL  = 364,BLKS IZE3368,DENx2) 

//SYS  IN  DD  * 

SENET  FRAME  OCCUPANCY 
5000  15000  1000 
TQUEUE  10000 
TVFRAME  15440 
TSFRAME  15440 
TNVOICE  193 
/* 


Figure  15.  JCL  For  Plot  Program 


similar  to  those  used  in  the  statistics  program  w|th  the  addition  of  a 
first  card  for  the  plot  title.  The  variables  are  normalized  before 
plotting.  Therefore,  the  maximum  value  input  scales  the  tabulated 
values  to  a realistic  range. 

The  JCL  for  translating  and  compiling  a model  is  given  in  Figure 
17.  The  source  is  stored  in  dataset  R3363. SOL. DATA  member  MODEL. 

The  translated  output  is  stored  in  dataset  R3363.S0L.PLI  in  a member 
of  the  same  name.  This  output  is  then  compiled  and  saved  in  dataset 
R3363.0BJ . Finally,  this  dataset  is  linked  and  the  load  module  is 
stored  in  R3363.L0AD(M0DEL). 


Figure  16.  Plotted  Results 


//R3363S0L  JOB  (AlOO, ,163,9,200) , 'CASTELLO‘ , 

//  NOT IFY=R3363 ,MSGCLAS$=Q 

//TRN  EXEC  P6M=CASTRAN , REG I0N=200K 

//STEPLIB  DO  0SN=UR1591.L0AD,DISP=SHR 

//SOLTRAN  DO  DSN=R3363.S0L.PL  KMODEL ) ,DISP*SHR 

//OUT FI  DD  SPACE=(CYL,(1 ,1) ) ,DISP= (NEW, DELETE) ,UNIT*SYSDA, 

//  DCB=(RECFM=FB,LRECL=80,BLKS IZE*3120) 

//SYSPRINT  DO  SYSOUT=(A,U) 

//TSOOUT  00  SYSOUT=(A,U) 

//SYS IN  DO  DSN=R3363. SOL. OATA( MODEL ),DISP=S HR 
//PL I EXEC  PGM=IELOAA,REG I0N=180K, 

//  PARM='M, OBJECT, NODECK ' ,C0ND=(EVEN,(0,LT,TRN) ) 

//SYSUT1  DD  UNIT=SYSDA,DCB=BLKS  IZE=1024,SPACE=(CYL ,(1 ,1 ) ) 
//SYSUT2  DD  UNIT=SYSDA,SPACE=(CYL ,(1,1)) 

//SYSLIN  DO  D ISP=S HR, DSN=R3363. OBJ (MODEL) 

//SYSPRINT  OD  SYSOUT=(A,U) 

//SYS  IN  DD  DISP=SHR,DSN=R3363.S0L.PL I(MOOEL) 

//SOL INC  DD  QISP=SHR,DSN=R1591 .CASTLIB.PL I 
//LKD  EXEC  PGM=IEWL,REGI0N=120K,C0ND=(EYEN) , 

//  PARM=,SIZE=(114688,26624) ' 

//SYSLIB  00  DSN-SYS1.PLIBASE,DISP=SHR 
//SYSLIN  OD  DSN=R3363.0BJ( MODEL) ,DISP=SHR 
//SYSLMOO  DO  0SN=R3363. LOAD( MODEL ) ,DISP=SHR 
//SYSPRINT  DD  SYSOUT=(A,U) 

//SYSUT1  DD  UNIT=SYSDA,XB=BLKS  IZE=1024,SPACE=(CYL ,(1,1)) 


Figure  17.  JCL  Translate,  Compile,  and  Link 
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APPENDIX  A 


SPEECH  ACTIVITY  MODEL  DEVELOPMENT 

This  appendix  discusses  the  development  of  the  speech  activity 
model  used  in  this  paper.  The  appendix  begins  with  a short  discussion 
of  the  speech  process  and  previous  authors'  experiments  and  models  of 

r 

speech.  Then  the  development  of  the  speech  model  used  in  this  report 
is  traced. 

■ 

Speech  is  not  an  orderly  process.  A speaker  utters  sound  for  a 
short  period  of  time  and  than  pauses  before  continuing  to  talk  again. 

Some  of  these  pauses  are  so  short  that  the  listener  does  not  notice 
the  discontinuity.  A "talkspurt"  is  speech  by  a person  which  is 
discerned  by  the  listener  as  continuous  speech.  Talkspurts  are 

terminated  when  the  speaker  drops  silent.  These  relatively  long  *1 

periods  of  silence  between  talkspurts  are  defined  as  "silent  periods." 

In  addition  to  silent  periods  and  talkspurts,  the  short  silences 
within  a talkspurt  are  of  interest.  These  short  silences  are  caused 
by  stop  consonants,  word  breaks,  and  other  interruptions  which  are 
called  "gaps"  in  this  report.  These  definitions  remain  the  same  for 
conversational  speech,  but  the  generating  mechanisms  change. 

Conversations  consist  of  interactive  speech  between  two  or  more 
people.  The  speech  is  interactive  in  that  speakers  interrupt  one 
another  and  alternate  talking.  Interruptions  are  either  short  monosyl- 
lables, such  as  "yes,"  "right,"  etc...  or  takeover  statements. 

Takeover  statements  force  the  intial  speaker  to  become  silent  while 


54 


. _ —*■  iitlf 


the  Interrupting  speaker  continues  talking.  Polite  conversation  also 
occurs  In  which  a second  speaker  waits  until  the  first  is  finished  and 
then  begins  to  speak.  These  interactions  are  the  difference  between 
monologue  speech  and  conversational  speech. 

Soaie  of  the  differences  between  conversational  speech  and  mono- 
logue speech  will  not  Impact  model  development;  others  will.  In 
conversation,  talkspurts  are  the  same  as  In  monologue  speech,  but  they 
may  terminate  because  of  intervention  by  the  other  speaker.  The  gaps 
within  a talkspurt  are  not  impacted  because  they  occur  due  to  the 
characteristics  of  the  language  and  the  physical  speech  process.  The 
silent  periods  will  be  modified  because  a speaker  will  now  pause  for 
the  same  reasons  as  in  monologue  speech,  plus  waiting  for  the  other 
speaker  to  pause.  This  change  to  silent  periods  also  influences  the 
arrival  rate  of  talkspurts.  Conversational  speech  also  Increases  the 
number  of  short  monosyllabic  replies.  We  need  to  keep  these  differ- 
ences In  mind  as  we  consider  speech  models. 

Measurements  have  been  made  of  both  monologue  and  conversational 
speech.  In  this  report  we  are  concerned  mainly  with  measurements  of 
conversational  speech.  The  best  documented  experiments  are  those  done 
by  Norwine  and  Murphy  [4]  and  Brady  [5].  While  other  authors  have 
hinted  at  experimental  work,  the  results  of  the  experimental  on  have 
not  been  published  directly. 

i £ 

Norwine  and  Murphy's  experiment  was  done  in  the  late  1930's.  Due 
to  the  type  of  equipment  available  to  them  and  the  lack  of  a detailed 
technical  description  of  the  speech  power  levels  involved,  a comparison 
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performed  their  experiments  to  enable  Horton  to  apply  probability 
theory  to  a study  of  the  occurrence  of  lockouts  caused  by  echo  suppres- 
sors and  similar  equipments  on  long  distance  telephone  circuits  [6]. 
Because  Brady's  measurements  indicate  that  the  statistics  of  speech 
events  change  with  detector  characteristics,  the  work  of  Norwine  and 
Murphy  is  only  of  historical  significance. 

Brady  published  a series  of  papers  in  the  Bell  System  Technical 
Journal;  the  first  two  described  his  experimental  setup  for  measuring 
speech  activity  [7]  and  the  results  of  measurements  of  16  conver- 
sations [5].  His  reason  for  studying  speech  activity  is  similar  to 
Norwine  and  Murphy's  - to  learn  more  about  conversational  speech 
patterns  prompted  by  the  introduction  of  long  telephone  lines  with 
significant  delay  and  the  increased  number  of  voice-actuated  devices. 
Brady's  equipment  detects  speech  by  rectifying  the  incoming  audio  and 
comparing  the  rectified  signal  to  a reference  level.  Each  time  the 
input  signal  exceeds  the  reference  level  a flip-flop  is  set.  A clock 
signal  resets  the  flip-flop  every  5 ms.  The  output  of  the  flip-flop 
is  then  recorded  on  magnetic  tape,  from  which  the  signal  is  then 
digitally  processed.  Note  that  with  this  equipment  any  signal  exceed- 
ing the  threshold,  no  matter  how  short  or  when  it  occurs  within  the  5 
ms  interval,  sets  the  flip-flop  for  the  entire  sample  duration. 

Digital  processing  of  the  magnetic  tape  created  by  this  sampling 
technique  removes  speech  detected  (active)  periods  less  than  or  equal 
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to  15  ms  and  speech  silence  periods  less  than  or  equal  to  200  ms. 
Therefore,  the  speech  experiment  does  not  consider  either  short  speech 
active  intervals  or  gaps  in  talkspurts. 

Using  this  procedure  Brady  discovered  that  speech  characteristics 
are  threshold  level  sensitive.  Some  of  the  results  of  his  measure- 
ments are  tabulated  in  Table  A-I.  Brady's  measurements,  however,  do 
not  give  us  insight  into  short  gaps  in  speech. 

TABLE  A-I.  MEASUREMENTS  OF  SPEECH  ACTIVITY 
Detector  level  in  dBrnU 

27.8  minutes 


-45 

-40 

-35 

test  period 

Talkspurts 

Percent 

43.53 

39.50 

35.00 

Mean 

1.366 

1.197 

0.98 

Number 

5486 

5794 

6224 

Mutual  silence 

(both  speakers 

silent) 

Percent 

18.97 

25.01 

32.55 

Mean 

1.802 

1.845 

1.742 

Number 

5485 

5792 

6424 

Sherman  in 

a concise  paper  to  the 

IEEE  Transactions  on 

Communica- 

tions  mentions 

work  he  performed  with 

Brady  which  measures 

the  gaps  in 

talkspurts  [8], 

A talkspurt 

is  broken 

i into  active  periods 

with  mean 

50  ms  duration  and  gaps  with  mean  10  ms. 

Stochastic  models  for  the  speech  process  developed  by  various 
authors  give  us  insight  into  the  development  of  a model  for  use  in 
this  paper.  Most  models  of  speech  that  this  author  has  found  are 
stochastic  models  in  continuous  time  for  a single  speaker  or  a single 
conversation.  Models  also  vary  in  the  number  of  states  used  to 
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represent  the  speech  process.  The  simplest  model  for  monologue  speech 
consists  of  two  states  and  the  most  complex  model  for  conversational 
speech  consists  of  six  states. 

Models  of  speech  activity  have  been  developed  by  many  authors. 

In  the  area  of  psychology,  Jaffe,  Cassotta,  and  Feldstein  [9]  devel- 
oped a simple  two  state  model  with  exponentially  distributed  hold 
times,  and  Gustafson  [10]  developed  a two  state  model  with  modified 
geometric  distributions  for  speech  activity.  In  the  case  of  a hundred 
simultaneous  voice  calls,  Birdsall,  Ristenbalt,  and  Weinstein  approxi- 
mated talkspurt  arrivals  with  a Poisson  model  [11].  Brady  [12] 
developed  a six-state  model  for  conversational  speech,  and  Sherman  [8] 
incorporated  gaps  within  talkspurts  in  his  model  for  monologue  speech. 

The  models  developed  by  the  last  two  authors  are  described  in  more 
detail  in  the  following  paragraphs. 

Brady's  six  state  model  with  is  the  most  ambitious  model  of  speech 
activity.  The  development  of  this  model  is  based  on  his  experimental 
measurements  described  above  in  Table  A-I.  Brady's  model  simulates  a 
conversation  between  two  talkers  from  talker  A's  viewpoint  and  is  based  on 
his  previous  experimental  measurements  [5].  If  the  model  is  split  in  half 
horizontally,  the  top  half  represents  A talking  and  the  bottom  half  repre- 
sents A silent.  Split  the  model  vertically  and  B is  talking  on  the 
right  half  and  silent  on  the  left.  Vertical  transitions  are  deter- 
mined by  talker  A,  and  horizontal  by  B.  Associated  with  each  transi- 
tion path  between  states  is  a parameter,  q-jj,  of  a Poisson  process. 
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A "pulse"  from  the  Poisson  process  terminates  the  state.  Stated  in 
other  terms,  q-jj  is  the  average  rate  in  times/second  that  a transi- 
tion to  state  i will  occur  given  that  the  system  is  in  state  j. 
Therefore,  speaker  A exists  in  the  one  of  the  six  states  and  his 
talk-silence  behavior  is  governed  by  a Poisson  process  with  time 
independent  parameters  and  by  the  state  he  occupies.  The  values  for 
the  exit  parameters,  q^j,  are  given  in  Table  A- II.  The  Markov 
transition  matrix  for  this  system  is  given  by  Table  A- III. 


TABLE  A- II 

:.  BRADY  MODEL  EXIT  PARAMETERS 

0.0 

2.07 

2.15 

2.24 

1.03 

0.0 

0.3 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.3 

0.79 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.79 

0.0 

2.15 

2.07 

1.03 

2.24 

0.0 

TABLE  A-III. 

MARKOV  TRANSITION  MATRIX 

FOR  BRADY  MODEL 

l_e-1.09t 

0.49e’4*22t 

0.51e-4*22t 

0.53e-3*27t 

0.47e-3.27t 

0.0 

0.275c-1 -09t 

l_e-4.22t 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

2_e-4 . 22t 

0.0 

0.0 

0.275c-1 

.09t 

0.725e-1*09t 

0.0 

0.0 

l_e-3.27t 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

l_e-3.27t 

0.725a-1 

• 09t 

0.0 

0.51e-4-22t 

0.49e-4.22t 

0.47e‘3*27t 

0.53e-3-27t 

l_e-1.09t 

Comparing  the  above  model  with  his  experimental  work,  Brady 
presented  the  following  conclusions: 

1.  Speech  activity  is  dependent  on  the  characteristics  of  the 
speech  detector. 

2.  Talkspurts  tend  to  have  an  exponential  distribution. 
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3.  Silent  periods  cannot  be  accurately  modeled  with  a simple 


ponentlal  distribution. 

solve  this  last  problem  Sherman  developed  an  alternative  technique 
r modeling  silent  periods. 

In  the  first  part  of  his  paper  [8]  Sherman  describes  an  alterna-  • 
ve  to  Brady's  modeling  of  silences.  He  noted  from  measurements  of 
nologue  speech  that  silent  periods  fall  into  two  general  catagories: 
art  silence  and  long  silence.  The  silent  periods  were  then  modeled 
a mixture  of  these  two  types. 

In  Sherman's  model  of  speech  activity  a talker  may  exist  in  one 
three  states:  talking,  long  pause  or  short  pause.  The  short 
jse  models  the  silent  periods  between  words  and  letters  within  a 
Ikspurt  and  the  long  pauses  represent  the  silent  periods  between 
itences.  This  model  allows  us  to  see  the  "fine"  structure  of 
lkspurts(See  next  page).  The  transition  matrix  for  this  model  may  be 
itten  in  terms  of  parameters,  q-jj,  as  in  the  Brady  model  above. 

Is  results  in  the  matrix  in  Table  A-IV,  and  the  steady  state  probabil- 
ies  shown  in  Table  A-V. 


TABLE  A-IV.  SHERMAN  MODEL  EXIT  PARAMETERS 

0.0  1.0  100.00 

1.0  0.0  0.0 

19.0  0.0  0.0 


— 1 

TALKING 


i 

> i 

’ 1 

2 

~*1 

LONG 

SHORT 

SILENCE 

SILENCE 

MONOLOGUE  SPEECH  MODEL 
by  SHERMAN 


SHERMAN  MODEL.  Figure  A- 2 


TABLE  A-V.  STEADY  STATE  PROBABILITIES 
FOR  SHERMAN  MODEL 


STATE 

1 

2 

3 


PROBABILITY 

.4566 

.4566 

.0868 


The  Markov  transition  matrix  for  this  system  is  given  by  Table  A-'/I. 
TABLE  A-VI.  MARKOV  TRANSITION  MATRIX  FOR  SHERMAN  MODEL 


l_e-20t 

e't 

e-100t 

0.05e-20t 

l-e"t 

0.0 

0.95e-2°t 

0.0 

l_e-100t 
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This  section  develops  a voice  activity  model  for  use  in  the 
analysis  of  SENET.  As  mentioned  previously,  SENET  handles  each 
direction  of  a full  duplex  circuit  independently.  Therefore,  only  one 
speaker  need  by  modelled.  The  monologue  speech  activity  model  of 
Sherman  modified  for  the  statistics  of  conversational  speech  well 
suits  our  requirements.  This  will  provide  us  with  a continuous  time 
model  for  speech  activity. 

The  conversational  speech  model  will  consist  on  three  states  ar- 
ranged in  the  same  form  as  in  Sherman's  model.  The  first  state  models 
the  actual  sounds  made  by  the  speaker.  It  is  assumed  that  the  duration 
in  this  state  may  be  modelled  by  an  exponential  density  with  mean  50  ms. 
The  second  state  represents  the  long  pauses  between  talking  states  due 
to  pauses  between  thoughts  and  pauses  while  the  other  speaker  is  talking. 
The  duration  in  this  state  Is  assumed  to  also  have  an  exponential  density 
with  mean  1/x.  The  third  state  which  represents  the  short  silences  within 
a talkspurt  remains  unchanged  from  Sherman's  model. 

The  most  significant  changes  from  Sherman's  model  is  in  the  duration 
of  the  silent  period  and  the  rates  at  which  states  two  and  three  are 
entered  from  state  one.  The  assumption  that  state  one  has  an  exponential 
duration  with  mean  50  ms  implies  that  the  exit  rate  from  this  state  will 
be  20  times/second.  If  we  let  y represent  the  rate  of  exit  to  state  two, 
then  state  one  must  exit  to  state  three  at  a rate  of  20-y.  From  the 
assumptions  made  for  states  two  and  three  they  will  have  exit  rates  to 
state  one  of  x and  100  respectively.  These  are  shown  on  Figure  A-3. 


Based  on  the  results  of  Brady's  experimental  measurements  values 
may  be  determined  for  x and  y.  Brady  has  Indicated  that  the  measure- 
ment made  at  -4S  dBm#  are  too  cluttered  with  noise  and  those  made  at 
-35  dBm#  clip  talkspurts.  Therefore,  we  will  use  the  statistics  of  voice 
measured  at  -40  dBm#.  Brady  filtered  out  the  very  short  talkspurts  and 
bridged  the  short  silences.  Therefore,  his  "talkspurt"  encompasses  what 
we  consider  as  two  separate  states:  states  one  and  three. 

If  we  consider  states  one  and  three  as  a single  talkspurt  state  then 
x is  the  rate  of  talkspurt  arrivals.  Based  on  Brady's  measurements  on 
conversation  speech  made  at  -40  dBm#  there  were  5794  talkspurts  during 
the  274.8  minutes  of  conversation  speech.  This  results  in  an  arrival 
rate  In  talkspurts  per  second  or  an  exit  rate  from  long  silent  periods  of 

x - 5794  * 0.3514. 

WTTTtt 

A value  for  y may  be  determined  from  Brady's  measured  probability 
of  being  in  the  silent  state.  Again  using  the  -40  dBm#  measurements 
the  probability  of  being  in  state  two  is  0.605.  We  assume  that  this  is 
the  steady  state  probability  of  being  In  our  state  two.  The  following 
method  may  be  used  to  determine  a value  for  y.  Let  p^  be  the  steady 
state  probability  of  being  in  state  1.  The  following  equations  hold 
for  steady  state. 
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„ -TZE  *5«5 ! 


yPi  = 0.3514  P2 

(20-y)p!  = 100p3 
1 = Pl  + p2  + P3 

P2  = 0.605 

In  matrix  notation  this  may  be  written  as 


y 

0.3514 

0 

PI 

0 

20-y 

0 

100 

0.605 

= 

0 

1 

1 

1 

P3 

0 

Solving  results  in  y = 19.37. 

With  these  values  for  the  unknowns  the  final  continuous  time 
model  is  given  in  Figure  A-4.  The  steady  state  probabilities  are 
given  in  table  A-V II,  the  transition  matrix  is  given  in  Table  A-VII 
and  the  continuous  model  exit  parameters  are  given  in  Table  A-IX. 

TABLE  A-VII.  STEADY  STATE  PROBABILITIES 
FOR  CONTINUOUS  MODEL 


STATE 


PROBABILITY 


SPEECH  MODEL.  Figure  A-4 


TABLE  A-VIII.  MARKOV  TRANSITION  MATRIX  FOR  CONTINUOUS  MODEL 
!_e-20t  e-0 .351 4t  e-100t 

0.0315e-20t  i_e-0.3514t  o.O 

_0.9625e"20t  0.0  l-e“100t 

TABLE  A-IX.  CONTINUOUS  MODEL  EXIT  PARAMETERS 


0.0 

0.3514 

100.00 

1.63 

0.0 

0.0 

19.37 

0.0 

0.0 
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